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ABSTRACT 

We estimated the luminosity function (LF) of IRAS galaxies in the PSCz catalogue. The faint end of the 
PSCz LF is slightly steeper than that of the LF derived by Saunders et al. (1990; S90). Using an analytical form 
for the LF used by S90, we obtain the following parameters: a = 1.23 ± 0.04, U = (8.85 ± 1.75) x 10 8 /T 2 L Q , 
a = 0.724 ±0.010, and <j)* = (2.34 ±0.30) x I0~ 2 h 3 Mpc" 3 . We also examined the evolution in the sample by a 
simple assumption <j>*(z) oc (1 +z) p , and found P = 3.40 ± 0.70. It does not affect the three parameters, a, L„, and 
a, but cj)*(z = 0) is overestimated up to ~ 15 % if we ignore evolution. We estimated the temperature dependence 
of the LF. The LFs of warm and cool galaxies are quite different: the LF of warm galaxies has a very steep faint 
end with a = 1 .37. We also discuss a lump found at the brightest end of the LF. 

Subject headings: galaxies: luminosity function, mass function — galaxies: statistics — infrared: galaxies — 
methods: statistical 



1. INTRODUCTION 

Galaxy luminosity function (LF) is an important tool to char- 
acterize the statistical properties of galaxies. The far-infrared 
(FIR) LF enables us to evaluate the properties of interstellar 
dust and its heating sources in various galaxy populations. The 
Infrared Astronomical Satellite (IRAS) has provided us a uni- 
form FIR all-sky survey of local galaxies. In addition to its own 
importance, the FIR LF has also come to the limelight in con- 
junction with the evolution of FIR galaxies. Evaluation of their 
evolution plays a crucial role in understanding the cosmic star 
formation history hidden by dust (e.g., Franceschini et al. 2001; 
Takeuchi et al. 2001a,b; Rowan-Robinson 2001; Granato et al. 
2000; Totani & Takeuchi 2002; Takeuchi, Shibai, & Ishii 2002). 

In spite of many attempts to estimate the LF of IRAS galax- 
ies (e.g., Lawrence et al. 1986; Soifer et al. 1987; Saunders et 
al. 1990 (S90); Isobe & Feigelson 1992), there is yet room for 
improvement in its exact shape, especially at the faint end. The 
best way to overcome this difficulty is, of course, a statistical 
analysis of a huge homogeneous dataset. 

The PSCz catalogue (Saunders et al. 2000) (S00) is the 
largest well controlled redshift sample of IRAS galaxies to date. 
The catalogue is complete to 0.6 Jy at 60 /im and contains 
1541 1 IRAS galaxies covering 84 % of the sky. This is an ideal 
dataset to estimate the exact shape of the LF of IRAS galaxies. 

In this Letter, we estimate the LF from the PSCz catalogue 
using the statistical package for LF estimation developed by us 
(Takeuchi, Yoshikawa, & Ishii 2000; hereafter TOO) and exam- 
ine the exact shape of the FIR LF. We use the cosmological pa- 
rameter set (h,Clo, Ao) = (0.7,0.3,0.7) except otherwise stated. 

2. ANALYSIS 

2.1. Data Description 

2.1.1. Flux density 

The data provided by S00 include (1) PSC flux, (2) point- 
source-filtered ADDSCAN flux, and (3) extended (coadded or 
extended ADDSCAN) flux for all the sources. First we should 



note that the flux densities of nearby extended sources can be 
seriously underestimated because of the fixed aperture of IRAS. 
Indeed, Figure 2 of S00 shows a clear trend of underestimation 
in PSC flux or point-source-filtered flux. We used (3), coadded 
or extended ADDSCAN flux, that has been carefully measured 
and homogenized by S00. 

2.1.2. Local Group galaxies 

The main catalog of S00 explicitly excludes the Local Group 
(LG) galaxies, and they are separately compiled in another cat- 
alog 'ilg.dat'. There is a possibility that the exclusion of these 
galaxies slightly affects the faintest-end slope estimate of the 
LF, because they generally have low luminosities. In order to 
include the LG galaxies we adopted the distances presented by 
van den Bergh (2000), adjusted for the cosmological parameters 
used. We compared the LFs estimated from the datasets with 
and without LG galaxies, and obtained completely the same re- 
sults: their inclusion caused no difference at all. For Virgo infall 
correction, we adopted a flow model very similar to model V3 
ofS90. 



2.2. K-correction 

To make ^-correction, we used extended fluxes of S25, S^o, 
and S100 tabulated in the PSCz catalogue and fitted a 2nd-order 
polynomial to the flux densities for each galaxy: 



log5 A (logA) = J>(logAy, 



logS™ = log%,o ( i +z)] 



(1) 



(2) 



where S\ is the interpolated flux density at wavelength A /im, 
superscript 'em' means that the flux is measured at the emitted 
wavelength, and coefficients a -, are estimated from the observed 
fluxes S25, $60, and Sioo- 
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2.3. Statistical Methods 

TOO extensively discussed four nonparametric statistical 
methods for estimating galaxy LF. The methods are 1)1 /V mm 
method (Schmidt 1968; Eales 1993), 2) Efstathiou, Ellis, & Pe- 
terson (1988) (EEP) method, 3) Choloniewski (1986) method, 
and 4) improved C~ method (Lynden-Bell 1971; Choloniewski 
1987). The associated error is estimated by a bootstrap resam- 
pling (see TOO). Equation (3) has been verified to be appropriate 
for a FIR LF (see Appendix D of S90). Hence, in Section 3 we 
also show the parametric fit for the LF by the parameter estima- 
tion advocated by Sandage, Tammann, & Yahil (1979), using 



(3) 
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We note that, while the C~, EEP, and the Choloniewski meth- 
ods are robust against density fluctuation, the 1 /V max method is 
not accurate if a density inhomogeneity exists in the sample, as 
pointed out by previous studies. 



significantly deviates upward because of the local density en- 
hancement. We note that LG galaxies are included in the sam- 
ple to produce Figure 1 . 

S90 proposed the following values for Equation (3): a = 
1.09 ± 0.120, U = (2.95! 2 ;^) x 10 8 /T 2 L 0) a = 0.724 ± 0.031, 
and 0* = (2.6 ±0.8) x I0~ 2 h 3 Mpc" 3 . These values were ob- 
tained from the compiled redshifts of 2818 IRAS galaxies. We 
find that the bright end of the PSCz LF is remarkably well de- 
scribed by the analytic form estimated by S90. The normaliza- 
tion also nicely agrees with that of S90. At this stage we ignore 
the possible effect of galaxy evolution, which we will discuss 
in Section 5.1.2. 

However, the faint end deviates from the result of S90. The 
three inhomogeneity-insensitive estimates show consistently 
the same trend. Hence the steep faint end of the PSCz LF is 
not caused by density enhancement, and we conclude that it is 
an intrinsic statistical property of the local IRAS galaxies. 

We also estimated parameters for the analytic fit expressed 
by Equation (3) as follows: 



3. THE PSCZ LUMINOSITY FUNCTION 
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FIG. 1.— The luminosity function (LF) of IRAS PSCz galaxies. The 
squares, diamonds, triangles, and crosses represent the results obtained with 
1/Vmax, improved C, EEP, and Choloniewski methods, respectively. The 
faint-end estimates (10 6 -10 9 Lg) are rebinned to reduce statistical fluctuation. 
The dashed line describes the analytic LF derived in this work. The parametric 
result of Saunders et al. (1990) is also shown by the dotted line. 

Figure 1 is the FIR LF of the PSCz catalogue. Different sym- 
bols represent the results obtained by different estimators. The 
dashed line describes the analytic LF derived in this work. The 
parametric solution of S90 is also shown in Figure 1 by the dot- 
ted line. We see that the estimates by the three inhomogeneity- 
insensitive methods show excellent agreement with one an- 
other, and at a faint regime Leo % 10 9 L Q , 1/V ma x estimate 



a =1.23 ±0.04, 

U =(8.85 ± 1.75) x 10 8 /T 2 [L ], 

a =0.724 ±0.01, 

0. = (2.60 ± 0.30) x 10" 2 /z 3 [Mpc" 3 ] . 



The uncertainty for each parameter is estimated from the rel- 
ative marginalized logarithmic likelihood. The log-likelihood 
around the maximum likelihood solution asymptotically be- 
haves like a Gaussian, hence we can estimate 1-er error from 
Aln£ = ln£-ln£ max = -0.5. We note that a is larger than that 
of S90 (3-sigma significant), though within their quoted error, 
it is consistent with the estimate of S90. The large a quan- 
titatively describes the above-mentioned steep faint end. On 
the other hand, we obtain almost the same value for a, which 
characterizes the bright-end shape, and <f>, which determines the 
density normalization. 



4. THE LUMINOSITY FUNCTION OF WARM AND COOL IRAS 
GALAXIES 

In this section we derive the LFs of relatively warm and cool 
galaxies in the PSCz sample of galaxies. We chose the bound- 
ary of these classes to be (3 = Siqq/S^q = 2.1, which is almost 
the median value for the sample. Helou (1986) has discussed 
the color-color diagram of normal IRAS galaxies. From his 
Figure 1, we find that (3 = 2.1 is also equal to the median of his 
sample. If we use a modified blackbody with emissivity index 
1.5, (i = 2.1 corresponds to the dust temperature <~ 25 K. We 
refer to the galaxies with (i < 2. 1 as 'warm galaxies', and those 
with f3 > 2.1 as 'cool galaxies' in this study. 

S00 have assigned extended ADDSCAN fluxes to all galaxies, 
we simply used these values to estimate /3. However, warmer 
galaxies tend to have larger luminosities (e.g., Takeuchi et al. 
2001a), this procedure might introduce a complex selection 
bias. We reconstruct the complete subsample of each category 
by adopting a shallower flux limit of S^o > 0.912 Jy. The final 
sizes of the two samples are 5554 (warm) and 3421 (cool). 
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FIG. 2. — The luminosity function (LF) of warm (J3 < 2.1) and cool 
(J3 > 2.1) galaxies. The solid curve describes the analytic expression of the 
LF of the whole PSCz galaxies. Filled and open squares represent the LFs of 
warm and cool galaxies, respectively. Dotted and dot-dashed lines depict the 
analytic fits of the warm and cool galaxies. 

The LFs of the warm and cool PSCz galaxies are shown in 
Figure 2. In Figure 2, we show only the estimates by EEP 
method for clarity; but the other estimates show a good agree- 
ment with one anther. The parameters of the analytic fit for 
warm galaxies are 



a = 1.37 ±0.05, 

U = (5.10 ± 0.90) x 10 9 /T 2 [Lq] , 

a =0.625 ±0.015, 

<f>* = (4.20 ± 0.85) x 10" 3 /z 3 [Mpc" 3 ] , 



(5) 



and those for cool galaxies are 



a = 1.25 ±0.06, 

L* = (1.95 ± 0.40) x 10 8 /T 2 [Lq] , 

a =0.500 ±0.020, 

4>* = (1.85 ± 0.37) x 10" 2 n 3 [Mpc" 3 ] . 



(6) 



We also show the analytic fits for these two groups in Figure 2. 
S90 performed a similar analysis of warm and cool galaxies. 
They first fit modified blackbody spectra to IRAS flux data and 
defined dust temperature, and then divide the sample into galax- 
ies with T > 36 K and T < 36 K. Though they give system- 
atically higher temperature than ours, their LF of each class 
is quite similar to ours; therefore, S90 and we probably treat 
nearly the same populations. 

It is interesting to note that the faintest end of the warm 
galaxies increases steeply. This means that a significant fraction 
of galaxies at luminosities Z^o ^ 10 7 Lq can be low-luminosity 
but hot galaxies. Because of their low luminosity, detailed case 
studies of their spectral energy distributions and next genera- 
tion large FIR surveys (e.g., ASTRO-F, SIRTF, and Herschel) 
are important to reveal their physical properties. 



5. DISCUSSION 

5.1. Possible Systematic Uncertainties in LF Parameters 

5.1.1. Hubble constant 

We further consider some subtle effects on the parameter es- 
timation. The effect of the choice of cosmological parameters 
is examined in S90. Their Figure 5 illustrates the dependence 
of the estimated LF. We should be cautious that the faint end 
slope is slightly affected by the choice of the Hubble constant, 
h. But according to their Figure 5, if we use larger h, the esti- 
mated a will be smaller. Hence, our steep faint-end slope is not 
due to this effect. 

5.1.2. Evolution 

Rowan-Robinson (2001) mentions the possibility that the pa- 
rameters are affected by galaxy evolution of the sample. Hence 
we examined this issue by dividing the whole sample into 
z < 0.02 and z > 0.02 groups. We derived the LF of each group 
and found that the faint end of the LF of the total PSCz sample 
is dominated by galaxies at z < 0.02. The flux density limit of 
0.6 Jy corresponds to L 60 ~ 10 9 ' 6 Lq at z = 0.02 (h = 0.7), there- 
fore the evolution toward higher redshift can hardly affect the 
estimation of a. 

On the other hand, evolution may affect the estimation of the 
normalization of the LF. We estimated the evolution strength 
under the assumption «(z) oc «o(l +z) p , where no and n(z) is the 
comoving number density of galaxies at redshift and z, re- 
spectively (S90). If we assume that the LF is separable for L 
and z as <j)(L,z) = n(z)p(L), then the likelihood is expressed as 



£(P\{(L i , Zi )} i= i,., N ) = ll 



n(zi)p(Li) 



i{ / n(z)p(Li)(dV/dz)dz 

TT (l+Zif 

\ = \j^'(l+z) p (dV/dz)dz 



, (7) 
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where p(L) denotes the probability density for a galaxy having 
a luminosity L, and Zmax,; represents the maximum redshift to 
which a galaxy ; can be detected within the flux limit Su m . We 
can estimate the strength of density evolution, P, by maximiz- 
ing the above likelihood equation [Equation (7)]. 

We obtained P = 3.40 ± 0.70. This value is much smaller 
than that estimated by S90, but consistent with the estimate 
of Springel & White (1998). By using this value, we found 
that the density normalization 0*(z = 0) is overestimated up to 
~ 10-15 % if evolution is ignored. This systematic error is 
comparable to the statistical uncertainty of ~ 12 %. Thus, the 
density normalization parameter is 0* = 2.34 if we properly take 
into account the evolution. 

5.2. Bright End Population 

It is worthwhile to note the lump found in the bright end of 
the LF in Figure 1 at Lgo ^ 10 L . The zoom-up of the 
bright-end is shown in Figure 3. We only plot EEP estimates 
for clarity. We consider which population of galaxies dominates 



the high-luminosity lump. Machalski & Godlowski (2000) pre- 
sented the radio (1.4 GHz) LF from UGC (Nilson 1973) and 
LCRS (Shectman et al. 1996) galaxy samples. Their radio 
LF clearly shows that it consists of two distinct populations: 
star-forming galaxies and AGNs. At Li.4gh z > 10 23 ' 4 [WHz -1 ] 
AGN contribution dominates the LF. We can convert the ra- 
dio luminosity to FIR one by the well-known radio-FIR cor- 
relation (Condon 1992) and obtain the corresponding L F tr ~ 
5.0 x 10 11 Lq. It agrees with the 60-/im luminosity at which 
the lump begins to appear. So we suggest that the AGN con- 
tribution is attributed to the lump of the LF at L^q > 10 11,5 L Q . 
Different energetics may be the cause of discontinuity in the LF 
shape. It is interesting that both warm and cool galaxies show 
departure from the analytic fit at the same L^q (see Figure 2). 

We thank two referees for thoughtful comments and scrutiny 
that have improved the quality of this paper very much. We 
also thank H. Hirashita, T. N. Rengarajan, and M. Imanishi who 
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